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Abstract 

We have searched for the lepton-flavor-violating decay t ^ using a data sample of 86.7 fb~^ 
collected with the Belle detector at the KEKB asymmetric e"'"e~ collider. No evidence for a signal 
is obtained, and we set an upper limit for the branching fraction B{t — > 67) < 3.9 x 10~^ at the 
90% C.L. 



PACS numbers: 13.35.Dx, 11.30.Fs, 14.60.Fg 



INTRODUCTION 



Lepton-flavor-violating (LFV) processes are good probes of physics beyond the Standard 
Model (SM). For instance, in some supersymmetric models, off-diagonal components of the 
left-handed slepton mass matrix, mi, could radiatively induce LFV such as in r — > /"(e)7 
and /X — >• 67 decays [1, 2]. In general, the branching fraction B{t — >• fij) is expected to be 
larger than B{t 67), since the mixing between the third and second families is typically 
assumed to be stronger than that between the third and first families. However, if the first 
and third families couple more strongly, for instance, due to an inverted hierarchy of slepton 
masses, then B{t 67) could exceed B{t — > n'y) and might be detectable [3]. Values of 
B{t — >■ 67) which can exceed that for B{t — > /X7) are also predicted in the models with 
heavy Dirac neutrinos [4, 5]. Thus, a study of both r — 67 and r — 1^7 decays is essential 
not only to search for new physics but also to further examine lepton flavor structure. 

The decay r — > 67 has been searched for, along with r — > /i7, by MARK II [6] , Crystal 
Ball [7], ARGUS [8], DELPHI [9], and CLEO [10], among which CLEO has set the most 
sensitive upper limit of B{t 67) < 2.7 x lO"*" at 90% C.L. 

Recently the Belle collaboration performed a search for the LFV decay r — > /X7 [11]. Here 
we present a new search for the decay r ^ 67 based on data samples of 77.7 fb~^ and 9.0 
fb~^, collected at the T(4S) resonance and in the continuum 60 MeV below the resonance, 
respectively, equivalent in total to 77.3 x 10® r+r~ pairs. The data were collected with the 
Belle detector at the KEKB asymmetric e^e~ collider [12]. A description of the detector 
can be found in Ref. [13]. 

DATA SELECTION 

We search for events containing exactly two oppositely-charged tracks and at least one 
photon. The events should be consistent with a r+r" event in which one r (signal side) 
decays to ej and the other (tag side) decays to a non-electron charged particle (denoted 
hereafter as ^), neutrino(s) and any number of photons. 

The selection criteria are determined by studying Monte Carlo (MC) simulations for 
signal r-pair decay and background (BG) events, such as generic r-pair decay (r+r"), qq 
continuum, BB, Bhabha, fj^'^^~, and two-photon events [11]. The KORALB/TAUOLA [14] 
and QQ [15] generators are used for event generation, and the Belle detector response is 
simulated by a GEANT3 [16] based program. The two-body decay r — > 67 is initially 
assumed to have a uniform angular distribution in the r lepton's rest system. 

The selection criteria are similar to those used in the r — /i7 search [11]. Kinematic 
variables with a CM superscript are calculated in the center-of-mass frame; all other variables 
are calculated in the laboratory frame. Before electron identification, all the charged tracks 
are assumed to be pions. Each track is required to have momentum < 4.5 GeV/c 
and momentum transverse to the e"*" beam pt > 0.1 GeV/c, the former requirement being 
imposed to avoid Bhabha and fi'^fi~ contamination. We require that the energy of 
each photon exceed 0.1 GeV. In addition, we also require the total energy measured in the 
CsI(Tl) electromagnetic calorimeter (ECL), E-^cl: to be less than 9 GeV in order to suppress 
background from Bhabha events. 

The tracks and photons must be detected within the detector's fiducial volume —0.866 < 
cos^ < 0.956, but outside the barrel-endcap gaps defined by 0.829 < cos 6* < 0.880 and 
—0.716 < cos^ < —0.602. Here, 9 is the polar angle with respect to the direction opposite 
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FIG. 1: (a) Energy distribution of the signal candidate photon, (b) cos^^^ distribution, (c) 
distributions. The open histogram is the sum of backgrounds from generic r+r~, qq (uds) contin- 
uum, Bhabha, and two-photon processes evaluated from MC simulation. Dots indicate the 
data distribution, and the shaded histogram is the signal MC distribution. Electron identification 
requirements were applied for these figures. 



to the e"*" beam. Identification of electrons is performed using an electron likelihood ratio, 
£e, which is based on the dE/dx information from the central drift chamber (CDC), the 
ratio of the energy deposited in the ECL to the momentum measured by both the CDC and 
the silicon vertex detector (SVD), the shower shape in the ECL, the hit information from 
the aerogel Cherenkov counter, and time-of-flight measurements [17]. 

The electron track that forms a r — > 67 candidate (hereafter denoted as (67)) is required 
to have an e likelihood ratio > 0.90 and a momentum p > 1.0 GeV/c. This requirement 
has an efficiency of (93±3)% in the barrel and forward detector and (76±7)% in the backward 
detector because of the additional material. On the tag side, the ^ track is required to have 
Ce < 0.1. The fraction r} of electrons with >Ce < 0.1 is measured to be (4 ± 3)% in the barrel 
and forward detector and (13 ±5)% in the backward detector for p > 1.0 GeV/c. 

The photon that forms an (67) candidate is required to have > 0.5 GeV in order to 
reduce spurious combinations of a low-energy 7 with an electron, see Fig. 1(a). 

A requirement on the cosine of the opening angle between the e and 7 of the (57) candi- 
date, 0.4 < cos6'g^' < 0.8, is particularly powerful in rejecting the generic T+r~ BG events 
(see Fig. 1(b)). The events in Fig. 1(b) that peak at cos^^^^ ~ 1, arise from electrons that 
radiate a photon when they interact in the SVD or in materials around it. The requirement 
-^sum < 9.0 GeV is imposed to reject Bhabha and production, where E^^^ is defined 

as the sum of the energies of the two charged tracks and the photon composing the (67), 
see Fig. 1(c). The opening angle between the two tracks in the laboratory frame is required 
to be greater than 90°. 

We define Pmiss ^is the residual momentum vector calculated by subtracting the vector 
sum of all visible momenta (of both tracks and photons) from the vector sum of the beam 
momenta. Constraints on the momentum and cosine of the polar angle of the missing 
particle are imposed: pmiss > 0.4 GeV/c and —0.866 < cos^miss < 0.956. To remove t+t~ 
BG events, we apply a requirement on the opening angle between the tagging track and the 
missing particle of 0.4 < cos6'^^g_. < 0.99. 




FIG. 2: (a) cos^^^g.^^^g distribution. For ^-tagged events, the distributions of signal MC (his- 
togram), generic t~^t~ MC (boxes), and ^(e7)data sample (open circles) are shown. For e-tagged 
events, the distribution for e(e7) data (closed circles) are also plotted. All requirements except 
the one for cos ^m^g-tag applied, (b) cos 9^ and (c) cos Oe distributions. These are polar angle 
distributions of the tag side track for ^{ej) and e{ej) data, respectively, where the solid histogram 
is for the events with cos ^m^s_tag > 0-99 and the dotted histogram is for 0.4 < cos ^miss-tag < ^-^^ 
(tagged by e or ^). 

The upper bound on cos9^^_^ is introduced to reject radiative Bhabha events in which 
one of the electrons forms an (67) candidate with a radiated photon and the electron on 
the tag side is misidentified as the ^ due to the electron identification inefficiency. By 
analyzing a Bhabha data sample, a large portion of such events is found to have a very 
small opening angle, cos6^^_^ ^ 1, and a polar angle peaking strongly forward, cos^^^ > 
0.8. Figure 2(a) shows the cos ^^^j^gg_^^g distributions wdth tag given by ^ or e for ^(67) or 
e{ej) modes, respectively, in the actual Bhabha data samples, and the signal and generic 
T+T~ MC data. Figures 2 (b) and (c) present the cos^tag distribution for ^(67) and e{e^) 
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FIG. 3: Event distribution in the rn^i^^ — Pmiss plane for (a) signal and (b) generic r+r MCs. 
The events within the two lines are accepted in the analysis. 

Bhabha data samples, respectively. The requirement, cos^^^^g.^ag < 0.99, reduces ^{e^) and 
e(e7) candidates that originate from radiative Bhabhas by 73% and 45%, respectively, while 
only slightly affecting the signal (97%) and generic r+r^ (99%) events. 

Finally, a condition is imposed on the relation between Pmiss and the mass-squared of a 
missing particle, m^^^^. The latter is defined as E^^^^ — p^jsg, where -Bmiss is 11-5 GeV (the 
sum of the beam energies) minus the sum of all visible energy and is calculated assuming 
the electron (pion) mass for the charged track on the signal (tag) side. We require Pmiss > 
-5 (cVGeV) X mj2„i33 - 1 (GeV/c) and p^iss > 1-5 (cVGeV) x m^j^g - 1 (GeV/c), where p^iss 
is in GeV/c and mmiss is in GeV/c^ (see Fig. 3). With this cut, 98% of the generic r+r" 
and 97% of the e+e~7 backgrounds are removed, while 69% of the signal events remain. In 
addition, most of the remaining BB, continuum, and two-photon events are rejected by this 
requirement. 

After these selection requirements, 224 events remain in the data, about 3 times fewer 
than in the r — > case. Since the inefficiency of electron identification is much smaller 
than that of the muon, the Bhabha BG is strongly suppressed in spite of its much larger 
cross-section than that of e~^e~ /i"'"/i~7. The r ^ e'-f detection efficiency is evaluated 
by MC to be 7.29%, about 40% smaller than that of r — > /ij, mostly because of the E-ecl 
requirement. 

True signal events will have an invariant mass (Minv) close to the r lepton mass and an 
energy close to the beam energy in the CM frame, i.e., AE = E^^ — i?bcam — 0- When 
deciding on our selection criteria, we excluded the signal region 1.68 GeV/c^ < Mi^v < 
1.85 GeV/c^ so as not to bias our choice of criteria (a "blind" analysis). Only after all 
requirements were finalized and the number of BG events estimated did we include this 
region and count the number of signal events. 
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FIG. 4: (a) AE' distributions for the data (dots) and the expected BG (curve and open histogram) 
in the bhnded region. The distribution for signal MC is the shaded curve. See the text for more 
detail, (b) Minv vs. AE distributions for the data (dots) and signal MC (shaded boxes). The ±5(7 
region is indicated by the dashed rectangle. 



RESULTS 

Background Evaluation 

To analyze the BG distributions, we define a region, named "grand signal region": —0.5 
GeV< AE < 0.5 GeV and 1.5 GeV/c^ < M^^^ < 2.0 GeV/c^, containing 90% of signal MC 
events passing all previous requirements. 

The resolution in AE and Minv is evaluated by MC: an asymmetric Gaussian reproduces 
the dominant part of the signal MC distribution with a^^J"^'^^ = (84.8 ± 1.2)/(36.0 ± 0.9) 
MeV and a^f''^^' = (25.7 ± 0.3)/(14.3 ± 0.2) MeV/c^, where ^^""/high ^eans the standard 
deviation at the lower/higher side of the peak. The peak positions are —6.2 ± 1.0 MeV and 
1776.0 ± 0.2 MeV/c^ for AE and Mi„^, respectively. 

A dominant source of BG is the process e'^e~ — > t'^t~^, in which the photon is radiated 
from the initial state: an (67) candidate is formed by the electron from the r — > euV decay 



and the initial state radiation photon, while the tag side r decays generically via a one- 
prong mode but not to an electron. From a 174 fb~^ sample of MC t^t~7 events we find 
j\jTT^ _ gg g _j_ events in the "grand signal region." 

The contribution from the process ^67 was described above and is evaluated as N^'^'^ ~ 
K X N'^^'^, where k = ri/{l — rj). Prom the data N'^'^'^ is found to be 68.0 ± 8.2 events and 
K is estimated to be 0.06 ± 0.03 from both the Bhabha data and MC samples taking into 
account the momentum dependence based on the momentum distribution of the signal MC 
events. Thus, we have 7V^'^''' = 4.3±2.0 events. 

Prom the MC simulation, no other process is expected to contribute to the background. 
Therefore, the expected BG in the "grand signal region" is 65.1 ± 5.9 events. 

The Miiiv and Ai? shapes of both types of BG events are empirically reproduced by a 
combination of Landau and Gaussian functions. 

For r"'"r~7. 
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Here a ^ AE - AEll^y^{Mi^^) and f3 = AE - A^^p^^kl^mv), where AE^peak denotes the peak 
position in terms of c x Minv + d for individual BG's. The parameters a, b, c, d, vi/}^ and cuyh 
are determined from MC for r"'"r~7 and from data for the e+e~7. 

The BG distribution can then be represented by the sum of the two BG components 
above as 

7VBG(Min„ AE) = 7V^^^(Minv, AE) + k x N''^{M,^,, AE). (3) 

Figure 4(a) compares the AE distribution in the 1.70 GeV/c^ < Mjnv < 1.82 GeV/c^ 
(±3(TMinv) I'egion for BG events expected from Eq. (3) (the solid curve) and the events 
obtained by interpolating the data distribution from both sidebands, 1.53 GeV/c^ < Minv < 
1.68 GeV/c^ and 1.85 GeV/c^ < Mi^y < 2.0 GeV/c^ (the open histogram). Good agreement 
between them is observed. 



Upper Limit for 15{t 67) 



After opening the blinded region, we find the AE and M^^^ vs. AE distributions that 
are shown in Figs. 4 (a) and (b), respectively. The number of surviving data events in 



the "grand signal region" is 60, in good agreement with the expected BG contribution of 
65.1 ± 5.9 events. 

In order to extract the number of signal events from the surviving sample, we apply an 
unbinned extended maximum likelihood fit with the likelihood function defined as 

J^^^;TrIli^3i + bBi), (4) 

i=l 

where N is the number of observed events, s and b are free parameters representing the 
numbers of signal and BG events to be extracted, respectively, and Si = S(M^^^, AE'^'^^) 
and Bi = B{M^^^^, AE^^'') are the signal and BG probability density functions for the i-th 
event. The function B{M-^^^^, AE) is taken from Eq. (3) while S'(Minv, A£^) is obtained by 
generating 10^ signal MC events. 

We apply this fit for s and 6 to a ±5(j region in Minv and AE: 1.65 GeV/c^ < 
Min^ <1.85 GeV/c^ and -0.43 GeV < AE <0.17 GeV. There are a total of 20 events 
in this region while 25.7 ± 0.3 events are expected from Eq. (3), and, when s is constrained 
to be non- negative, the fit finds s = and h = 20.0. 

To calculate the upper limit, Monte Carlo samples are generated by fixing the expected 
number of BG events (b) to the value b — 20. For every assumed expected number of signal 
events {§), 10,000 samples are generated, for each of which the numbers of signal and BG 
events are determined by Poisson statistics with means § and b, respectively. We then assign 
Minv and AE values to these events according to their density distributions. An unbinned 
maximum likelihood fit is performed for every sample to extract the number of signal events 
^gMC^_ The confidence level for an assumed s is defined as the fraction of the samples whose 
gMC exceeds s. This procedure is repeated until we find the value of s (sgo) that gives a 90% 
chance of being larger than s. 

The resulting upper limit at 90% C.L. is Sgo = 3.75 events. An upper limit on the 
branching fraction is obtained via the formula: 

where A^,-^ is the total number of r-pairs produced, and e is the detection efficiency in the ±5a 
region. Inserting the values Nt-t- = 77.3 x 10^ and e = 6.37% gives B{t — > ej) < 3.8 x 10~^. 



Systematic Uncertainties 

Systematic uncertainties on Sgo are evaluated by varying all parameters of the BG proba- 
bility density function. The fractions of N'^'^"' [Mi-^^^ AE) and N'^'^'^{Mi^i^^ AE) in Eq. (3) are 
varied by ±20% and ±100%, respectively, about double their estimated uncertainties. As a 
result, Sgo varies by +0.01/— 0.00 and +0.01/— 0.02 events, respectively. The functional form 
of the BG spectra is scaled by 1.15 or 0.90 times for N'^'^'^ [Mi-ay, ^E) and by 1.3 or 0.6 times 
for A^^^T(Minv, AE), and their centers are shifted by +0.01/-0.015 GeV for Ar^^T(Minv, AE) 
and by ±0.1 GeV for N^'^^{M\a-y, AE), all changes corresponding to the estimated errors of 
the involved parameters. The shift of the central value for the N'^'^"' [M\-av, AE) spectrum 
yields the largest effect of +0.07/— 0.13 events, and the overall systematic uncertainty in- 
creasing the upper hmit of Sgo is evaluated as +0.07 events. The stability of the result for the 



fit region is examined by extending the Minv-A£' region from ±4(T to iScr: no appreciable 
difference in the upper hmit is found. 

The systematic uncertainties on the detection sensitivity, 2eA^^^, arise from the photon 
reconstruction efficiency (3.0%), the selection criteria (2.5%), the trigger efficiency (2.0%), 
the track reconstruction efficiency (2.0%), the luminosity (1.4%), and the MC statistics 
(0.3%). The total uncertainty is obtained by adding all of these components in quadrature; 
the result is 5.0%. The contribution of the largest component, the photon reconstruction 
efficiency, is evaluated from the e+e^7 data sample. The uncertainty of the selection criteria 
is estimated by varying the required polar angle region of the signal candidate photon. The 
trigger efficiency is estimated from the difference between a t^t~ data sample and a generic 
T^T~ MC sample. 

These uncertainties are included in the upper limit on B{t 67) following [18]. The 
systematic uncertainty in the efficiency is assumed to have a Gaussian distribution. 

While the angular distribution of the r — > 67 decay is assumed to be uniform in this 
analysis, it is sensitive to the LFV interaction structure [19], and spin correlations between 
the r leptons on the signal and tag sides must be considered. To evaluate the maximum 
possible variation, V — A and V + A interactions are assumed; no statistically significant 
difference in the Minv-A£' distribution or in the efficiency is found compared to the case 
of the uniform distribution. Therefore, systematic uncertainties due to these effects are 
neglected in the upper limit evaluation. 

The incorporation of all systematic uncertainties increases the upper limit by 2.1%. As 
a result, the upper limit on the branching fraction is 

B{t e7) < 3.9 X 10"^ at 90% C.L. (6) 

SUMMARY 

This result improves the sensitivity to the branching fraction by approximately one order 
of magnitude compared to previous measurements. Despite a smaller detection efficiency 
compared to r — /i7, the superior BG rejection for electrons allows us to reach a sensitivity 
for T — > 67 that is comparable to r — > /X7. 
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